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CHAOTIC  ELECTRON  MOTION  CAUSED  BY  SIDEBANDS  IN 
FREE  ELECTRON  LASERS 

I.  INTRODUCTION 


Multifrequency  effects  in  Free  Electron  Lasers  (FELs)  become 
increasingly  important  as  progress  is  made  towards  high  power  operation. 
Growth  of  parasitic  frequencies  (sidebands1-5)  has  been  predicted 
theoretically  and  has  been  observed  in  experiments®’ ^  as  well  as  in 
simulations®-1*  with  either  constant  or  tapered  wigglers11.  The 
efficiency  for  the  carrier  signal  is  reduced  and  the  optical  quality  is 
degraded  as  power  is  channeled  into  frequencies  apart  from  the  intended 
operation  frequency.  Another  potential  hazard  that  has  attracted  little 
attention  so  far  is  the  onset  of  chaotic  electron  motion  caused  by  the 
presence  of  even  a  single  frequency  sideband.  This  may  lead  to  extensive 
particle  detrapping  and  premature  loss  of  the  amplification  for  all  the 
radiation  modes  independent  of  frequency. 

Two  of  the  main  issues  concerning  FEL  operation  are:  (a)  whether 
unstable  parasitic  frequencies  exist  that  can  grow  to  significant 
amplitude  and  (b)  what  is  the  effect  of  potentially  unstable  modes  on 
the  trapped  electron  trajectories.  Considerable  attention  has  been 
devoted  to  the  linear  stability  issue.  The  gain  for  small  sideband 
signal  has  been  computed  analytically1-5  invoking  either  ensemble 
averaging  over  single  particle  trajectories  or  solutions  of  the 
perturbed  kinetic  equation  for  the  distribution  function.  Initial 
results,  obtained  for  particles  localized  near  the  bottom  of  the 
ponderomotive  well,  and,  in  particular,  more  recent . results  including 
all  trapped  and  untrapped  particles  with  arbitrary  distributions,  have 
demonstrated  that  every  nontrivial  distribution  dfQ/dJ  ^  0  is  unstable 
to  sideband  growth. 
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Given  that  sidebands  cannot  be  eliminated,  the  growth  of  the 
unstable  modes  to  a  finite  amplitude  may  have  serious  effects  on  the 

12 

unperturbed  trajectories.  It  has  been  known  that  stochastic  behavior 

13  14 

is  an  intrinsic  property  of  perturbed  Hamiltonian  systems  ’ 

Accordingly,  the  electron  motion  in  a  FEL  will  become  chaotic  when  the 

sideband  amplitude  exceeds  a  certain  threshold.  This,  in  turn,  will 

result  in  significant  electron  detrapping.  Since  it  is  the  deceleration 

of  the  trapped  electron  bucket  that  provides  the  energy  for  the 

radiation  in  case  of  tapered  wigglers,  detrapping  will  cause  loss  of 

amplification  for  the  FEL  signal. 

In  the  present  work  we  investigate  the  nonlinear  effects  caused 

by  sidebands.  The  threshold  for  stochasticity,  above  which  unbound 

chaotic  motion  occurs,  is  determined.  Once  the  stochastic  transition 

takes  place,  the  action  J,  a  constant  of  motion  in  the  unperturbed 

2 

system,  changes  in  a  random  manner.  The  ensemble  average  <AX  >  of  any 

physical  quantity  X  is  described  by  a  diffusion  equation.  Diffusion  of 

the  action  invariant  provides  a  measure  of  the  leakage  rate  across  the 

separatrix.  If  D  is  the  effective  diffusion  coefficient  in  action  space 

2 

then  the  diffusion  length  L ,  =  J  /D,  where  J  is  the  action  at  the 

d  s  s 

separatrix,  signifies  the  length  over  which  approximately  half  of  the 
deeply  trapped  particles  get  detrapped.  We  show  that  a  single  frequency 
sideband  at  a  modest  fraction  of  the  carrier  amplitude  suffices  to 
spread  irregular  motion  over  a  significant  fraction  of  the  trapped 
particle  domain.  However,  given  that  the  interaction  time  of  an  electron 
in  a  FEL  is  short,  we  are  mainly  concerned  on  how  fast  this  diffusion 
occurs.  The  diffusion  rate  increases  and  the  diffusion  length  L^ 
decreases  with  increasing  sideband  ampli tude(s) .  Thus,  a  critical 
sideband  level  ac  can  be  defined  above  which  the  diffusion  length  L^ 
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becomes  shorter  than  the  wiggler  length  Ly.  Obviously  the  power  level 
for  the  sidebands  in  a  FEL  cannot  exceed  ac,  otherwise,  extensive 
diffusion  and  premature  detrapping  will  occur.  On  the  other  hand,  enough 
electrons  can  remain  trapped  during  the  interaction  period  even  though 
their  motion  has  turned  stochastic,  because  we  find  that  usually  the 
critical  amplitude  ac  is  much  larger  than  the  threshold  for 
stochasticity  ag. 

A  clear-cut  relation  between  the  diffusion  rate  under  constant 
total  sideband  power  and  the  type  of  the  excited  sideband  spectrum  is 
discovered.  More  specifically  we  observe  three  regimes  in  the  simulation 
parameters  defining  the  sideband  spectrum,  corresponding  to  a  narrow,  a 
wide  discrete  and  a  wide  continuous  spectrum.  The  transition  from  one 
spectral  type  to  another  is  accompanied  by  an  abrupt  change  in  the 
diffusion  rates.  In  all  cases  we  find  the  diffusion  coefficient 
proportional  to  the  ratio  of  the  total  power  in  the  sidebands  to  the  FEL 
carrier  power.  The  coefficients  of  this  proportionality  depend  on  the 
spectral  type.  A  general  conclusion  is  that  the  diffusion  rate  under 
constant  sideband  power  ratio  decreases  with  increasing  number  of 
spectral  components.  The  diffusion  rate  for  a  single  sideband  frequency 
exceeds  that  of  a  broad  continuous  spectrum  by  orders  of  magnitude  while 
a  broad  discrete  spectrum  causes  intermediate  diffusion  rates. 

For  practical  purposes  we  measure  the  diffusion  length  in  terms 
of  the  number  of  wiggler  periods,  while  <AJ‘>  is  normalized 

to  the  action  J  at  the  (unperturbed)  separatiix.  We  compute  the 
norma  1 1  zed  diffusion  coefficient  D  analytically  lot  t ho  cases  of  broad 
discrete  and  continuous  spectra.  In  the  latter  case  the  quasilinear 


diffusion  coetilcient  in  action  space  D^(J)  is  obtained  in  closed  form. 
This  expression  for  D^(J)  is  quite  general,  valid  for  any  choice  of 


unperturbed  Hamiltonian  Hq(J).  The  analysis  also  shows  that  the 
normalized  diffusion  coefficient  does  not  depend  on  the  beam  energy  y^,. 
The  numerical  results  agree  well  with  the  theory. 

We  evaluate  the  loss  of  trapped  particles  for  typical  short 
wavelength  FEL  parameters.  Ue  find  that  a  single  frequency  sideband  with 
a  sideband  to  carrier  power  ratio  of  <  1  can  cause  half  of  the  particles 
to  detrap  over  100  wiggler  periods;  we  have  observed  total  loss  of 
trapped  particles  for  power  ratios  of  =1.  In  cases  of  wide  but  discrete 
sideband  spectrum  the  diffusion  length  becomes  comparable  to  the  wiggler 
length  only  at  large  power  ratios  (>  1).  The  case  of  a  wide  continuous 
spectrum  seems  to  cause  insignificant  electron  detrapping  for  the  same 
parameters  as  above;  the  typical  diffusion  length  is  of  the  order  of 
1000  wiggler  periods  for  sideband  to  carrier  power  ratios  of  1. 

In  our  investigation  we  have  assumed  all  electromagnetic  fields 
as  given.  The  changes  in  the  particle  trajectories  are  decoupled  from 
the  evolution  of  the  fields.  At  the  expense  of  self-consistency  we  are 
able  to  analyze  the  situation  theoretically  and  determine  the  scaling  of 
the  diffusion  rates  on  the  various  FEL  parameters.  Deterioration  in  the 
extraction  efficiency  has  been  observed  in  self-consistent  numerical 
simulations  of  high  power  FEL  oscillators^  with  high  level  sideband 
excitation.  The  gain  per  pass  in  a  tapered  wiggler  is  progressively 
limited  as  the  sideband  power  goes  up  and  the  rate  of  electron 
detrapping  is  accelerated.  In  an  untapered  wiggler,  on  the  other  hand, 
particle  detrapping  is  not  so  important  for  the  main  signal  efficiency. 
The  total  extraction  efficiency  may  actually  increase  with  the  sidebands 
since  there  are  more  modes  to  channel  the  electron  beam  energy  into. 


The  remainder  of  this  paper  is  organized  as  follows.  In  Sec.  II 
we  construct  our  analytic  model  for  the  study  of  the  stochastic 
diffusion  and  discuss  the  various  approximations.  To  elucidate  the 
analysis  we  start  with  a  single  sideband  mode  and  give  a  sketchy 
description  of  how  this  can  lead  to  electron  detrapping.  In  Sec.  Ill  we 
examine  the  structure  of  the  phase  space  for  a  monochromatic  sideband  in 
detail,  using  canonical  formalism.  The  threshold  for  the  stochastic 
transition  and  the  extent  of  the  chaotic  regime  in  phase  space  are 
obtained  in  Sec.  IV.  In  Sec.  V  the  diffusion  rate  caused  by  a  single 
sideband  mode  is  examined  in  connection  with  the  various  FEL  parameters. 
In  Sec.  VI  the  study  is  extended  to  broad  (multifrequency)  sideband 
spectra.  A  distinction  is  drawn  between  continuous  and  discrete  spectra. 
Subsection  VI. a  covers  the  case  of  a  broad  discrete  spectrum  and  the 
related  diffusion  coefficient.  Subsection  VI. b  deals  with  a  broad 
continuous  spectrum  and  the  corresponding  quasilinear  diffusion 
coefficient.  In  Sec.  VII  the  theoretical  models  are  compared  with 
numerical  results.  The  differences  in  the  induced  diffusion  rates  among 
the  three  different  types  of  spectra  are  emphasized.  The  reduction  in 
the  extraction  efficiency  in  a  tapered  wiggler  FEL  is  computed  as  a 
function  of  the  diffusion  coefficient.  Results  and  conclusions  are 


summarized  in  Sec.  VIII. 


II.  GENERAL  CONSIDERATIONS 


We  consider  relativistic  electrons  streaming  along  the  z- 
direction  through  the  static  magnetic  wiggler  and  the  radiation  fields 
of  the  carrier  and  the  sideband.  We  take  all  fields  to  be  circularly 
polarized  and  of  constant  amplitude.  To  simplify  the  analysis  and  make 
the  underlying  ideas  clearer  we  start  out  with  monochromatic  waves  for 
the  carrier  and  the  sideband.  The  total  vector  potential  is  then, 


A  (z,t)  = 
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where  the  subscripts  w,  r,  and  s  stand  for  wiggler,  carrier  and  sideband 
respectively.  We  assume  that  all  waves  propagate  with  the  speed  of  light 
c,  ignoring  the  small  correction  of  order  /<*>r  from  the  dielectric 
contribution  of  the  beam.  Electrostatic  contributions  to  the  fields  are 
neglected  for  operation  in  the  Compton  regime.  The  phase  of  the  wiggler 
is  given  by  <f»v(z)  =  J  kv(z')dz',  where  the  wave  number  kv( z )  may  change 
slowly  in  z  on  a  scale  length  much  longer  than  the  wiggler  wavelength  Ay 
=  2n/kv*  The  main  signal  wave  number  kr  is  doubly  Doppler  upshifted  from 
the  wiggler  wavenumber  ky, 


k  =  2y  2k  ,  (2) 

r  z  w  v  ' 

2  -1/2 

with  the  upshifting  factor  y  =  (1  -  6  )  and  6  =  «  /c(k  +k  ). 

z  r  r  r  r  v 

We  have  ignored  variations  in  the  x-  and  y-directions.  Increased 
number  of  dimensions  is  known  to  facilitate  the  transition  to  chaotic 


6 


motion.  r'.erefore,  the  threshold  for  stochasticity  for  variations  in  the 
z-direction  only  will  be  useful  in  providing  a  neccessary  condition  to 
avoid  fast  large-scale  diffusion.  Slow  diffusion  due  to  higher 
dimensionality  will  in  fact  persist  for  the  real  system  below  this 
threshold.  As  far  as  particle  detrapping  is  concerned,  three  dimensional 
effects  are  comparatively  insignificant,  provided  that  the  dependence  on 
x  and  y  is  adiabatic.  This  requires  that  the  frequency  of  the  betatron 
oscillation,  caused  by  the  transverse  field  gradients^,  be  small 
compared  to  the  electron  synchrotron  frequency  in  the  ponderomotive 
bucket . 

We  have  also  assumed  that  the  radiation  amplitudes  remain 
constant.  In  case  of  fast  growth  rate  of  the  carrier  amplitude  the 
particle  trajectories  are  not  analytically  tractable,  even  in  the 
absence  of  sidebands.  It  is  generally  expected  that  the  fraction  of 
trapped  particles  decreases  with  decreasing  carrier  amplitude.  Therefore 
the  spreading  of  the  radiation  beam  due  to  diffraction^’  ^  can  also 
cause  detrapping  by  reducing  the  carrier  amplitude  a^.  This  detrapping 
mechanism  is  independent  of  the  diffusive  detrapping  caused  by  sidebands 
that  is  examined  here. 

Normalizing  the  time  t  to  wr”*»  the  length  z  to  kr_*,  the  mass  to 

2 

mg  and  the  vector  potentials  according  to  a^=|e|A^/mec  the 
dimensionless  Hamiltonian  describing  the  electron  motion  in  the  fields 
of  Eq.  (1)  is, 


[2  2 

u  +  p  -  2A  a  a  cos($  +  kz-wt)  -  2A  a  a  cos(<f>  +  k  z  -  w  t) 

H  *z  wr  wr  r  ws  ws  s  j 


1/2 


2  .  u/  2  2  2 . 
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(3) 


with  M  =  1  and  A  =  1.  Eq.  (3)  also  describes  the  fast-time  averaged 
Hamiltonian  for  a  linearly  polarized  wiggler  by  setting  M  =  1/2, 

A  =  Uo(C)  -  J^Ol/2  and  C  =  ay2/(4  +  2ay2). 

The  terms  proportional  to  ayar  and  ayas  are  the  ponderomotive 
potentials  due  to  the  combined  action  of  the  viggler  with  the  main 
signal  and  the  sideband  respectively.  The  resonant  velocities  for  each 
ponderomotive  potential  are  given  by  3^  =  k^/(k..  +  ky)  corresponding  to 
resonant  energies, 


In  the  vicinity  of  y^  the  motion  of  the  electrons  is  determined 
by  the  corresponding  resonant  term  inside  (3).  We  may  drop  the 
nonresonant  term  for  small  radiation  amplitudes  and  linearize  (3)  for 
small  excursions  &y  around  y^,.  From  the  resulting  pendulum  equation  we 
find  that  trapped  electrons  will  undergo  oscillations  of  frequency 
around  y^,  forming  islands  of  width  Sy^  in  phase  space,  where  and  8y^ 
are  given  by, 


We  call  these  islands,  due  to  the  direct  wave-particle  resonances, 
primary  islands. 

Roughly  speaking,  irregular  motion  breaks  out  as  a  result  of 

14 

nearby  island  overlapping  .  The  amplitudes  a^,  ag  must  increase  to  the 
point  where, 


*Yr  ♦  *TS  >  |  Tr  -  rs|  , 


(6) 


for  an  overlapping  between  the  two  primary  islands  to  take  place.  The 

2 

difference  in  resonant  energies  Ay  =  y  -  y  is  given  by  Ay  =  y  y  A3 

r  s  z  it 

where  Ag  =  |g  -  g  j  =  (l/2y  2) jk  -k  j/k  .  Given  that  typically  k  -k  - 
sir  z  r  s  r  sr 

2 

2yz  (a^  we  find, 

Ay  ~  yr(awar)1/2.  (7) 

It  then  follows  from  (5)  to  (7)  that  overlapping  and  transition  to 

chaotic  behavior  can  take  place  at  a  -  a  .  This  crude  estimate 

s  r 

demonstrates  the  potential  of  chaotic  behavior  for  large  amplitude 
sidebands.  The  above  threshold  becomes  even  smaller  in  case  of  a 
multifrequency  sideband  spectrum. 

We  will  be  interested  in  evaluating  the  fraction  of  the  phase 
space  that  becomes  chaotic  as  a  function  of  the  sideband  amplitude.  This 
requires  the  use  of  a  more  refined  overlapping  criterion.  Electrons 
trapped  inside  the  primary  island  of  the  main  signal  still  experience 
perturbations  in  their  motion  caused  by  the  sideband.  The  perturbation 
is  especially  felt  by  these  electrons  that  have  the  synchrotron 
frequency  matching  the  difference  between  the  frequencies  of  the  main 
signal  and  the  sideband.  This  condition  defines  new  secondary  resonances 
between  the  electrons  and  the  sideband.  It  is  the  overlapping  among  the 
nearby  secondary  islands,  formed  inside  the  primary  island,  that 
determines  more  accurately  the  break  out  and  extent  of  the  stochastic 
behavior. 
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III.  PORTRAIT  OF  THE  PHASE  SPACE 


The  electrons  are  injected  into  a  FEL  with  energies  near  the 
resonant  energy  rr  for  the  main  signal  Expression  (3)  can  be 
linearized  for  small  excursions  Sy/Yr  <<  1  for  electrons  not  too  far 
from  the  separatrix.  Introducing  r  =  r  -  Yr  and  ♦  =  (ky  +  kr)z  -  w^t  as 
a  new  pair  of  canonical  variables  and  approximating  the  time 
t(z)  =  z/c0r  we  obtain, 


kw  ~2  awar  awas 

=  —  Y  +  -3 — (cosg/  +  ♦  singi  )  +  — cos(g/  -  5  z). 

Yr  r  r  r  Y  r  s 


(8) 


In  (8)  the  phase  flow  is  parametrized  by  the  traveled  length  z  inside 
the  wiggler  rather  than  the  time  t.  It  was  also  assumed  that  the  wiggler 
parameters  change  slowly  compared  to  the  wiggler  wavelength  2n/kv.  The 
term  sim|<r  parametrizes  the  rate  of  change  for  the  resonant  energy 
caused  by  the  change  in  the  wiggler  wavelength, 


d 

dz^r 


k  a  a 
r  w  r 


sin*r  , 


(9) 


where  g/  =  it  corresponds  to  an  untapered  wiggler.  The  term  5  in  the 
r  s 

sideband  phase  is  the  Doppler  downshifted  difference  between  the  signal 
and  the  sideband  wave  numbers, 


& 


s 


-  k  ) 
r 


(10) 
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In  the  absence  of  sidebands,  a  =0,  the  Hamiltonian  H  is 

s  o 

integrable.  The  unperturbed  trajectories  in  the  ponderomotive  well  are 
given  by  =  K.  where  K  is  the  reduced  energy  in  the  ponderomotive 

frame.  These  orbits  take  the  simplest  possible  form  expressed  in  terms 
of  the  action-angle  variables  (J,  9),  defined  as, 


J  =  ^  <fd*  Y(K,¥),  9  =  |jJ  d*'y(K,*')  ,  (11) 

where  K  =  Hq(J)  and  the  path  of  integration  is  over  the  unperturbed 

orbits.  For  trapped  particles  in  closed  trajectories,  the  action  J  is 

related  to  the  area  in  phase  space  enclosed  by  the  orbit.  For  untrapped 

particles  in  open  trajectories,  the  path  of  integration  in  Eq.  (11) 

depends  on  the  wiggler  type.  In  case  of  an  untapered  wiggler,  the  orbits 

are  periodic  and  the  limits  of  V  integration  are  from  0  to  2n.  In  case 

of  a  tapered  wiggler  the  path  of  integration  is  the  segment  of  the 

trajectory  that  begins  and  ends  at  enclosing  the  separatrix. 

Thus,  J  remains  finite,  avoiding  an  infinite  jump  in  action  across  the 

separatrix  that  would  result  by  considering  the  full  orbit  length  for 
18 

unbound  orbits  .  J  is  always  periodic  in  ty,  J(yt'I'1)  =  J(y»^2^  ^or  = 
^2  +  2n,  even  when  Hq  is  not  (case  of  tapered  wiggler). 

Hamiltonian  (8)  is  now  transformed  under  the  canonical 
transformation  defined  by  Eq.  (11)  into, 


a  a  00 

H( J , 9; z)  =  H  (J)  +  -2-£  E  Qn(J>  cos(n0  +  5  z )  +  Q  (J)  cos(n0  -  &  z). 

o  Yr  n=0  n  s  n  s 

(12) 


Qn+(J)  are  the  Fourier  coefficients  obtained  by  the  decomposition  of  the 
perturbing  sideband  phase  q/( J , 0)  -  5gz  into  harmonics  of  the  angle  9, 
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where  <KJ,  8)  is  obtained  by  inverting  Eq.  (11).  In  case  of  constant 
parameter  wiggler  J,  9  and  Qn  (J)  are  expressed  in  closed  forms  given  in 
Appendix  A. 

Hq(J)  is  independent  of  0  so  the  unperturbed  orbits  in  (J,  0) 
space  are  straight  lines, 

J  =  const.,  9=9  +  K.  ( J ) z  . 

O  D 

The  synchrotron  wave  number  «b(J)  is  connected  to  the  bounce  length 
and  the  synchrotron  frequency  in  the  laboratory  frame  (^(J)  with  the 
relation, 

dH  (J)  2n  «.(J) 

K.  (J)  =  — -  =  -  =  -2 -  .  (13) 

dJ  L  b(J)  c(3z 

Since  c  =  1  in  the  normalized  units  and  0z  =  1  in  the  cases  of  interest, 
we  may  use  (^(J)  in  place  of  «b(J)  as  well. 

Expression  (12)  for  the  transformed  Hamiltonian  reveals  the  new 
resonances  emerging  when  a  sideband  is  turned  on.  Defining  the  phase  of 
the  nth  sideband  induced  harmonic  9^  =  n0  ±  6gz,  the  stationary  phase 
condition  reads, 

±  nxb(J)  -  Sg  =  0,  or  +n0z<«>b(J)  -  =  0.  (14) 

Thus,  particles,  originally  in  unperturbed  orbits  J  =  Jn,  resonate  with 
the  sideband  when  the  nth  harmonic  of  their  synchrotron  period  wj)('ln) 
matches  the  downshifted  frequency  difference  between  the  sideband  and 
the  carrier  signal. 


For  a  given  n  and  sufficiently  small  a^  we  may  keep  only  the 

resonant  term  9^n^  to  examine  the  motion  in  the  vicinity  of  J  .  This  is 

n 

formally  achieved  by  the  canonical  transformation, 


9  =  n9  -  &  z  , 
s 


Z  =  Ssz  ,  1  =  K  +  i  J  , 

s 

coming  from  the  generating  function  F(0,z,I,Iz)  =  (n0  -  z)I  -  zl 
resulting  Hamiltonian  is, 


(15) 


Z*  The 


a  a 

Hn(I,Iz,e,Z)  =  Ho(nI)  +  ig(Iz  -  I)  +  Qn(nl)  cos6  +  0(asZ)  . 

(16) 

The  fixed  points  (Jn,  ©n)  are  found  from, 


d0 

dZ 


3H 


31 


3H 
_ o 

31 


(17a) 


dl 

dZ  " 


3H 
_ n 

39 


a  a 
w  s 


Qn(nl)  sin9  =  0. 


(17b) 


Using  relations  (15)  for  the  transformed  variables  we  recover  from  (17a) 
the  resonant  condition  (14)  while  (17b)  indicates  9n  =  kit  /  n,  k  = 
0,1,..., n— 1 . 

In  short,  a  single  frequency  sideband  causes  chains  of  secondary 
islands  to  appear  inside  the  original  primary  island.  Each  chain 
corresponds  to  a  given  harmonic  n  and  is  centered  around  the  stable 
fixed  points  J  ,  0  .  The  structure  of  the  phase  space  is  shown  in  Figs. 

1  and  2.  They  are  surfaces  of  section,  created  by  numerically 
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integrating  the  original  equations  of  motion  from  Hamiltonian  Eq.  (8) 
and  then  recording  the  intersection  point  of  each  trajectory  with  the 
plane  z  =  2fi/Sg.  The  y  vs.  ^  plots  are  on  the  left  side  in  Figs.  1  and 
2.  The  plots  on  the  right  side  show  the  same  surfaces  of  section  in 
action-angle  variables,  produced  by  the  transformations  (11).  The  bounce 
frequency  around  a  given  secondary  island  is  found  by  linearly  expanding 
the  resonant  Hamiltonian  (16)  in  61  =  I  -  In«  From  the  resulting 
pendulum  equation  and  from  relations  (15)  one  finds  that  the  secondary 
synchrotron  period  Qn  near  the  center  is  given  by, 


n 


while  the  half-width  of  the  island  8J  is 

n 


Wn  "  [2 


VV 

<dVdJ)J 


(18b) 


Representation  (12)  for  the  Hamiltonian  (8)  is  formally  independent  on 
the  details  of  the  transformations  (11).  Consequently,  the  same 
stability  analysis  applies  for  constant  as  well  as  variable  parameter 
wigglers. 
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IV.  THRESHOLD  FOR  ERRATIC  MOTION 


When  the  sideband  amplitude  exceeds  a  certain  amplitude  ag 

regarded  as  the  stochasticity  threshold,  the  presence  of  even  one 

sideband  frequency  suffices  to  transform  the  regular  coherent  motion, 

such  as  in  Fig.  1,  to  the  irregular  unbounded  motion  shown  in  Fig.  2. 

The  mechanism  for  this  radical  change  in  behavior  can  be  briefly 

described  as  follows.  The  trajectories  emanating  from  the  unstable  fixed 

points  (X-points)  of  a  secondary  island  do  not  actually  join  smoothly 

12  13 

around  that  island.  They  intersect  infinite  times  with  each  other  ’ 

due  to  the  effect  of  the  other  harmonics  n'  t  n  that  were  ignored  during 

the  local  approximation  Eq.  (16).  A  thin  layer  of  fuzzy  motion  thus 

surrounds  each  island  chain  of  given  n.  As  the  amplitude  ag  increases, 

the  width  of  each  island  increases  according  to  (18b)  and  so  does  the 

thickness  of  the  stochastic  layer  around  that  island.  At  a  given  point 

the  stochastic  layers  around  the  two  neighboring  island  chains  n  and  n+1 
14 

overlap  ,  allowing  particles  to  hop  from  one  island  to  another.  This 
signifies  the  beginning  of  unbounded,  random  motion  in  J  characterized 
as  stochastic  diffusion. 

Various  methods  of  different  accuracy  have  been  developed  for 

12  14 

estimating  the  stochasticity  threshold  ’  .  An  approximate  criterion 

that  works  well  in  most  cases  is, 

SJ  +  SJ  ,  >  |  AJn  ,  (19) 

n  n+1  -  3  n 

where  SJ  ,  SJ  .  are  the  separatrix  half-widths  and  AJ  =  J  ,  -  J  is 
n  n+1  n  n+1  n 

the  distance  between  the  separatrix  centers  for  the  n  and  n+1  harmonics 
respectively.  For  small  widths  SJ  and  distances  AJ  compared  to  J  we  may 
expand 
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t 


(20) 


,,  .  ,,  .  f3“b)  W 

“b(Jn.l>  -  W  "  I5rj  Mn  =  n 

n 

and  use  (18b)  with  J  ,  =  J  to  obtain  the  amplitude  a  (n)  for 
'  n+1  n  r  s 

overlapping 

.  .  (2  )2  W  (  Yr  )  . 

as(n)  "  UnJ  (do^/dJJj  hawQn(Jn)J  *  (21) 

n 

The  outermost  islands  centered  at  J  =  J  ,  correspond  to  larger  shear 

ns 

d«k(J)/dJ,  smaller  ^(J)  and  higher  harmonics  n,  for  given  a>s  -  w^. 

According  to  (21)  the  threshold  ag(n)  is  lower  near  the  separatrix  and 

the  outermost  secondary  islands  will  be  the  first  to  overlap.  The 

overlapping  is  progressively  extending  to  smaller  J  and  lower  n  as  a 

n  s 

increases.  The  macroscopic  stochastic  layer  first  appears  near  the 
original  separatrix  of  the  primary  island  and  spreads  to  the  interior  of 
the  trapped  particle  bucket.  We  take  the  amplitude  when  the  two 
innermost  harmonics  overlap  as  the  threshold  for  "global"  stochastic 
transition,  ag  s  ag(n^).  The  lowest  possible  harmonic  n^  for  given 
frequency  is  defined  by  the  resonant  condition  (14).  For  Jr  small  we 
have  both  &J  and  AJ  of  order  J  and  the  approximations  that  led  to  (21) 
are  not  valid.  In  this  case,  the  exact  expressions  for  J  and  w^(J)  must 
be  applied  inside  the  criterion  (19). 

We  may  obtain  the  dependence  of  ^(n)  on  the  various  parameters 
using  Eqs.  (A4)  and  (A7),  setting  d/dJ  =  (dX/dJ)  d/dX  and  utilizing  the 
properties  of  the  elliptic  integrals  to  compute  the  derivatives  in  Eq. 
(21).  We  find  that, 


i6 


1 


(22) 


as<n> 


2 

n 


F2(Xn) 


f 


where 


F(Xn) 


E,  (X  )  X  2  (X  2-  1) 
lv  ir  n  n  ' 


E2<Xn) 


(1  -  E1<V 


The  threshold  for  extensive  stochasticity  ag  a  as(n^)  is  independent  of 

y  and  a  .  The  trapping  parameter  X  is  determined  uniquely  from  J 
v  n  n 

according  to  X  =  H  (J  )y  /a  a  +  1/2  (see  Appendix  A).  Thus,  the 
n  o  n  r  v  r 

sideband  frequency  u>s,  related  to  ^(J^  through  the  resonant  condition 

(14),  is  the  only  parameter  that  «s/ar  depends  on.  The  scaling  in  (22) 

is  still  valid  in  case  of  large  secondary  island  width  with  a 

modification  in  the  numerical  factor  F. 

In  Fig.  3a  we  plot  in  solid  line  the  threshold  a  for  extensive 

stochasticity,  when  the  two  innermost  secondary  island  chains  overlap, 

as  a  function  of  the  frequency  difference  co^  -  &>r.  The  dotted  line  shows 

the  threshold  for  overlapping  between  the  next  two  secondary  island 

chains.  Some  deeply  trapped  orbits,  near  the  center  of  the  original 

primary  island,  still  persist  when  a  is  close  to  a  .  The  extent  of  the 

s  s 

area  unaffected  by  the  irregular  motion  when  ag  =  ag  is  given 
approximately  by  J  <  Jc  where  Jq  =  Jn  -  &Jn-  In  Fig.  3b  we  plot  the 
portion  Jc/Js  of  the  remaining  "good"  trajectories  when  the  sideband 
amplitude  equals  ag  as  a  function  of  <x>s  -  «r.  It  is'  seen  that  the 
threshold  as  is  larger  and  the  extent  of  the  stochastic  regime  is 
maximized  as  well  for  frequency  mismatch  near  a  harmonic  of  the 
synchrotron  frequency  ^(0)  at  the  bottom  of  the  ponderomotive  well. 


X  t 


The  twieshold  os  for  overlapping  is  considerably  lower  but  the  extent  of 
the  stochastic  regime  also  diminishes  for  frequencies  far  from  a 
harmonic  of  the  central  synchrotron  frequency. 

A  typical  phase  portrait  for  a  sideband  amplitude  ag  slightly 
above  ag  is  shown  in  Fig.  2b.  Two  different  kinds  of  regions  coexist:  a 
stochastic  regime  where  diffusive  behavior  prevails,  interrupted  here 
and  there  by  islands  of  regular  motion,  remnants  of  the  original  regular 
motion.  The  stochastic  regimes  are  interconnected  allowing  unbounded 
particle  transport.  The  rate  of  diffusion  as  well  as  the  decorrelation 
times  are  not  uniform  in  phase  space  but  depend  on  both  J  and  9. 

When  a  is  increased  well  above  a  the  chaotic  motion  engulfs 
s  s 

almost  100%  of  the  phase  space  (Fig.  2c).  The  decorrelation  time  is 
short  everywhere  in  phase  space.  In  this  parameter  regime  the  behavior 
of  the  system  can  be  described  by  a  diffusion  coefficient  D(J)  depending 
on  the  action  J  only  and  insensitive  to  the  frequency  of  the  driving 
sideband.  Total  stochastization  of  the  island  interior  occurs  roughly 
when  the  sideband  amplitude  grows  to  the  point  where  the  stable  fixed 
point  y  =  0,  <1/  =  It  +  at  the  center  of  the  original  island  becomes 
unstable. 


18 


V.  NARROW  FREQUENCY  BAND  DIFFUSION 


We  will  examine  first  the  diffusion  caused  by  the  presence  of  one 

single  frequency,  large  amplitude  sideband.  This  is  a  relevant 

approximation  in  case  of  a  narrow  sideband  spectrum.  The  term  narrow 

implies  a  spectral  width  D«s  much  smaller  than  the  frequency  separation 

2 

w  -  «  ,  typically  of  the  order  2y_  w,  .  We  examine  the  evolution  of  a 
s  r  zd 

monoenergetic  distribution  f(J;z=0)  =  5(J  -  Jq)  by  numerically 

2  2 

integrating  the  equations  of  motion.  We  plot  <AJ  >,  <J>  and  2  <AJ  >/z 

against  the  distance  z  in  Figs.  4a,  4b,  and  4c  respectively.  The 

electrons  are  initially  uniformly  distributed  in  9  with  constant  action 

Jq  =  0.7  Jg.  Different  curves  in  the  same  frame  correspond  to  different 

sideband  amplitudes  ag  at  a  given  frequency  a>s. 

For  a  constant  diffusion  coefficient  D,  independent  of  J,  the 

2  2  2 

average  <J>  and  the  mean  square  deviation  <AJ  >  =  <J  >  -  <J>  would 
2 

evolve  as  <AJ  >  =  (1/2)  D  z,  <J>  =  JQ.  The  dashed  curves  in  Figs.  4a  - 

4c  correspond  to  a  sideband  amplitude  ag  below  the  stochastici ty 

2 

threshold  ag.  The  deviation  <AJ  >  asymtotes  to  a  constant  after  an 

2 

initial  increase  while  the  ratio  <AJ  >/z  tends  to  zero  for  large  z.  In 
this  case  stochasticity  is  localized.  Different  stochastic  regimes  are 
still  separated  by  "good"  integrable  orbits  (KAM  surfaces)  located  in 
between.  Electrons  diffuse  until  they  are  stopped  at  the  boundaries  of 
the  stochastic  regimes  that  "compartmentalize"  the  phase  space.  The 
solid  curves  in  Figs.  4a  -  4c  correspond  to  sideband  amplitude  above  the 
global  stochasticity  threshold  ag.  This  means  that  the  last  good  orbit 

has  been  destroyed  allowing  different  stochastic  regimes  to 

2  2 
interconnect.  <AJ  >  now  increases  monotonically  and  the  ratio  <AJ  >/z 

remains  finite  for  large  z.  The  fact  that  the  diffusion  rate  is  not 
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constant,  and  that  the  average  <J>  changes  away  from  the  initial  value 
Jq,  shows  that  D  depends  strongly  on  J. 

In  principle,  one  could  determine  a  local  D(J)  by  advancing  test 

distributions  6(J-J  )  of  various  J  over  short,  distances  z.  Then  the 

o  o 

Fokker-Planck  equation  for  any  initial  distribution  fQ(J)  could  be 
solved  numerically  using  D(J).  Here,  instead,  ve  elect  to  measure 
directly  the  effective  diffusion  rate  associated  with  a  given  type  of 
initial  distribution.  We  do  so  by  integrating  numerically  the  equations 
of  motion,  Hamiltonian  (8),  for  a  number  of  particles  (typically  400) 
assuming  constant  amplitude  for  the  electromagnetic  fields.  A  uniform 
initial  distribution  in  phase  space  with  trapped  particles  inside  the 
(unperturbed)  separatrix  is  chosen,  fQ(J)  =  [1  -  S(J  -  Jg))/Js  where  S 
is  the  step  function.  This  situation  is  relevant  with  the  operation  of 
tapered  wiggler  FELs  where  the  trapped  particles  in  the  ponderomotive 
bucket  are  decelerating,  falling  quickly  behind  the  untrapped  particles 
and  thus  creating  large  distribution  gradients  near  the  separatrix. 

The  two  questions  of  practical  interest  are  (a)  what  percentage 

of  the  particles  will  eventually  get  detrapped  and  (b)  how  fast  do  they 

leak  outside  the  separatrix.  For  our  uniform  initial  distribution  the 

maximum  fraction  of  particles  becoming  detrapped  equals  the  fraction  of 

the  inside  of  the  separatrix  area  that  becomes  chaotic.  In  Fig.  5a  we 

plot  the  fraction  f^  of  the  particles  that  cross  the  original  separatrix 

J  as  a  function  of  the  traveled  wiggler  length  for  values  of  q  =  a  /a 
s  S  L 

below  the  threshold  for  extended  stochasticity .  In  all  cases  an  initial 
stage  of  quick  diffusion  is  followed  by  a  long  period  where  the  average 
number  of  untrapped  particles  remains  practically  constant.  The  results 
are  consistent  with  the  existence  of  a  boundary  in  phase  space  (KAM 
surface)  separating  two  regimes:  the  one  of  unbounded,  chaotic  motion 
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from  the  one  filled  with  regular,  coherent  orbits  of  particles  that 

remain  trapped.  Only  electrons  in  the  area  between  the  last  integrable 

surface  and  the  old  unperturbed  separatrix  will  diffuse  until  that  area 

is  depleted.  A  fraction  1  -  f^  of  the  original  primary  island  area  will 

remain  trapped  for  an  arbitrarily  long  time,  as  long  as  ag  remains  below 

the  threshold  ag  associated  with  the  particular  sideband  frequency.  This 

fraction  is  shrinking  as  ag  increases  and  the  bucket  "peels  off".  The 

situation  when  a  exceeds  a  is  shown  in  Fig.  5b.  The  fraction  of 
s  s 

untrapped  particles  f^  reaches  1  in  all  cases,  meaning  complete  absence 

of  particle  confinement  in  the  bucket.  All  particles  can  eventually 

escape  with  a  rate  that  increases  with  increasing  ag. 

Numerical  results  showing  the  fraction  of  detrapped  particles  f^ 

after  100  wiggler  periods  as  a  function  of  ag/ar  are  plotted  in  Fig.  6 

for  various  sideband  frequencies  |»  -  w  |.  The  length  over  which 

s  r 

approximately  half  of  the  initially  trapped  particles  get  detrapped  will 
be  discussed  in  the  next  section,  in  comparison  with  the  diffusion  rates 
from  other  types  of  sideband  spectra. 
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VI.  BROAD  FREQUENCY  BAND  DIFFUSION 


So  far  stochastic  electron  detrapping  caused  by  a  single 

frequency  sideband  has  been  examined.  It  was  argued  that  when  the 

2 

excited  sideband  spectrum  is  narrow  enough,  i.e.,  Dws  <<  2yz  w^,  the 

situation  can  be  reasonably  approximated  by  a  single  frequency  sideband. 

Here,  we  consider  the  situation  when  a  broad  spectrum  of  frequencies 

2 

have  been  excited,  D<og  >  2yz  We  will  make  a  distinction  between  a 

continuous  and  a  discrete  spectrum.  In  case  of  a  discrete  spectrum  the 
distance  between  two  nearby  sideband  frequencies  is  much  larger  than  the 
width  of  an  individual  spectral  line.  In  the  opposite  case,  when  various 
peaks  in  the  spectrum  merge  together,  we  will  talk  about  a  continuous 
spectrum.  We  may  model  numerically  both  cases  by  introducing  a 
modulation  in  the  sideband  phase  of  Hamiltonian  Eq.  (11), 


k  a  a  a  a 

w  2  v  r  .  w  s  . 

H(  Y,  ^;z)  =  — •  y  +  - (cos*  +  *  sin *  )  +  - co s(*  +  A  sin  vz 

Yr  Yr  r  Yr 


Ssz)  ’ 


(23) 


that  is  transformed  in  action-angle  variables  as 


H( J , 0; z )  -  Hq(J)  +  J , 6; z) , 


Hj( J, 0;z)  = 

cl  3  00  oo 

—  s  Y.  J  (A)  V  Q+(J)  cos[n0  +  S  (m)z)  +  Q”(J)  cos[n0  - 
y  m  _  n  s  n 

'r  m=-°°  n=0 


(24) 

6s(m)z  ] . 


The  frequency  mismatch  values  6s(m)  anc*  tbe  corresponding  sideband 
frequencies  ws(m)  are  given  by, 
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(25) 


2 

6  (m)  =  5  +  m  v,  «  (m)  =  w  +  2  m  y  \>, 
s'  '  o  s'  7  so  'z  ’ 

2 

where  «  =  w  +  2y  h  .  Since  the  Bessel  function  coefficients  become 

so  r  z  o 

vanishingly  small,  Jm(A)  <<  1  for  A  >  m,  the  width  of  the  spectrum  is 

2 

given  by  D8  -  A  v  or  Da  =  2y  Av. 

s  s  z 

In  order  to  examine  the  connection  between  diffusion  rates  and 
the  types  of  the  sideband  spectra,  we  divide  the  latter  into  three 
general  categories:  narrow,  broad  discrete  and  broad  continuous.  The 
passage  from  one  regime  to  the  other  is  not  gradual  but  characterized  by 
abrupt  changes  in  the  diffusion  coefficients.  Thus,  from  the  diffusion 
point  of  view,  the  distinction  among  the  spectral  types  is  not  arbitrary 
but  based  on  certain  relations  between  the  parameters  A  and  v.  In  all 
three  regimes  of  the  parameter  space  the  rate  of  diffusion  is 
proportional  to  the  ratio  of  the  total  sideband  power  to  the  carrier 
power.  The  scaling  of  the  coefficients  of  this  proportionality  on  the 
various  FEL  parameters,  however,  differs  from  one  regime  to  the  other. 

Both  cases  of  the  broad  spectrum  are  characterized  by  a  width  Du>s 
in  the  excited  frequencies  that  is  larger  than  the  upshifted  synchrotron 
frequency  o^, 


9  n 

Do>s  >  2  >  equivalent  to  v  >  ^—  ,  (26) 

with  A  >>  1.  The  further  distinction  between  discrete  or  continuous 
spectrum  is  related  to  the  separation  between  nearby  frequencies.  We 

>  \>  >  w^/A  the  diffusion  rate  agrees  well  with 

the  quasilinear  diffusion  coefficient.  A  different  coefficient  is 

1/2 

derived  for  the  case  when  v  >  w^/A  >  “^/A,  in  agreement  with  the 

1/2 

numerical  simulations.  Consequenity  the  separation  w^/A  between 

2  i 


find  that  when  w^/A 


nearby  modes  marks  the  transition  from  a  discrete  to  a  continuous  type 

of  behavior.  Departure  from  the  quasilinear  diffusion  coefficient  has 
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also  been  observed  numerically  in  previous  literature  when  the 
frequency  separation  between  nearby  modes  was  not  "small  enough".  Here  a 
condition  for  the  discrete-to-continuous  transition  has  beeen  obtained. 

For  any  spectral  type,  the  sideband  amplitude  must  be  above  the 
stochastisity  threshold  in  order  to  trigger  electron  diffusion.  Using 
the  same  method  of  nearby  resonance  overlapping  as  in  Sec.  Ill,  and  the 
Hamiltonian  Eq.  (25)  we  find  that  the  threshold  in  case  of  a 
multifrequency  spectrum  is  given  by, 


“s  *  Vm<A)  “s  ’ 


(27) 


where  a  is  the  threshold  for  the  single  sideband  frequency.  Thus  a 
s  s 

decreases  with  decreasing  frequency  separation  v.  Condition  (27) 

guarantees  the  stochastization  of  the  particle  orbits.  The  frequency 

separation  among  sidebands  must  be  limited  by  the  additional  condition  v 
1/2 

<  w^/A  ,  as  stated  in  the  previous  paragraph,  if  one  wishes  to 

simulate  quasilinear  diffusion  with  a  discrete  spectrum. 


A.  Broad  Discrete  Spectrum. 


We  now  evaluate  the  diffusion  coefficient  for  a  broad,  discrete 
spectrum.  The  equation  of  motion  for  J  can  be  written  as, 


dJ  9H.  8H1  a  a  ,  >.  .  ... 

—  =  -  -  =  - - =  A  -  sin  +  A  sinvz  -  6  z  |V(\|/  )-  V ( \p)  |  , 

dz  89  dip  d9  y  V  S  '  mX 
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using  Eq.  (A2)  for  d^/d0.  Due  to  the  presence  of  many  frequencies  in  the 

spectrum  J(z)  executes  a  complicated  oscillatory  motion  with  the  average 

<  J  >  changing  very  little  most  of  the  time.  J  however  receives  a  large 

kick  AJ  near  resonances,  where  the  phase  *  =  +  A  sin  vz  -  &  z  of  the 

s 

multifrequency  perturbation  varies  slowly.  The  resonant  condition  is, 


d* 

dz 


+  A  v  sin  vz  -  =  0, 


(28) 


at  some  z  =  z^.  Given  that  kwy/Yr  <  collective  effects  due  to  many 
frequencies  are  important  for  the  resonance  in  Eq.  (28)  when  A  v  > 

On  this  basis  inequality  (27)  signifies  the  transition  from  a  narrow  to 
a  broad  spectrum.  Let  us  consider  the  case  A  v  »  w^.  Then  the 
resonances  occur  at  z^  =  in/v,  i  integer,  and  the  interval  between 
successive  resonances  is  Az  a  n/v.  Expanding  the  phase  #(z)  in  the 
equation  of  motion  for  J  around  the  resonance  Zj, 


♦(z)  a  +  j  ^2(sini</i  +  sin^r)  +  A  v2cos  vz.j  (z  -  Zj)2,  (29) 


and  extending  the  limits  of  integration  to  z  =  ±®  we  obtain, 


DJi  = 


V>r  (  2^1/2  cos  (*„  ±  J)  |v<*mx)  -  V(*i)|1/2  _ 

2  1/2  ' 

Yr  “b  [  (sin^  +  sin^r)  +  A  (~“)2) 


where  ^mx(J)  is  the  turning  point  for  an  unperturbed  trajectory  of  given 
2  2 

J.  When  A  v  >>  both  Az  and  AJ  depend  on  the  features  v  and  A  of  the 

sideband  spectrum  and  not  on  the  bounce  frequency  We  classify  the 

1/2 

cases  with  frequency  separation  v  >  w^/A  as  broad  discrete  spectra. 
They  obey  a  distinct  scaling  in  the  diffusion  coefficient  that  will  be 
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1/2 

derived  below.  We  find  the  cases  with  v  <  co^/A  to  agree  numerically 
with  the  quasilinear  diffusion  that  will  be  studied  in  the  next 
subsection. 

The  resonant  phases  between  two  successive  jumps  of  AJ  become 
quickly  decorrelated  when  ag  grows  above  the  stochastisity  threshold. 
Because  the  relation  between  and  involves  the  distance  z^,  and 
ij>.  will  also  become  decorrelated,  <cos*^  cos^>  =  0.  Assuming  complete 
decorrelation  between  two  successive  jumps  we  obtain, 


2  <4J,>  2  a  2  a  2 

D.  -  — - >■  <31> 

Dz  y  A  v 

r 

where  the  angular  brackets  <...>  signify  the  average  over  for 

constant  J.  For  practical  purposes  it  is  more  convenient  to  rescale  the 

diffusion  coefficient  so  that  the  distance  z  =  z/X  is  measured  in  terms 

—  w 

of  wiggler  wavelengths  and  the  action  J  =  J/Jg  signifies  the  location 
relative  to  the  separatrix.  In  these  units,  using  Eq.  (A4)  for  Jg  and 
setting  v  /  w^O)  =  r  we  obtain, 


D 

v 


3  2, 

n  g  C 


8  A  r 


(  Vr 


n1/2 


U  +  a. 


(32) 


The  term  g  is  a  scaling  factor,  the  ratio  of  the  separatrix  area  for  an 
untapered  wiggler  to  that  of  a  tapered  wiggler,  g  =  Js(»j»r*rt)/J(V>  ),  and 
depends  only  on  The  term  C  ~  1,  coming  from  the  averaging  over  tfc 
in  Eq.  (31),  is  computed  in  Appendix  B.  The  typical  diffusion  length  L^, 
the  traveled  distance  inside  the  wiggler  over  which  the  average  trapped 
particle  crosses  the  separatrix,  is  estimated  from  the  diffusion  Eq. 
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I 

(22)  by  talcing  <M2>  =  J  2  L,  ~  J  2/D  .  Thus,  the  diffusion  length  in 

S  Q  S  W 

wiggler  periods  =  L^/Xy  is  the  inverse  of  D^, 


B.  Broad  Continuous  Spectrum 


Next  we  consider  the  case  of  a  sideband  wave  package, 


as(z,t) 


1_ 

2n 


I 


dk  a  (k  ) 
s  s  s 


ik  z  -  i«(k  )t 
s  s 


(34) 


of  finite  spectral  width  Dk  centered  around  k  .  Our  purpose  is  to 

s  so 

obtain  the  diffusion  coefficient  for  a  continuous  spectrum  using  the 
methods  of  the  quasilinear  theory.  Upon  using  exression  (34)  for  the 
fields,  the  Hamiltonian  representation  in  action-angle  variables  assumes 
the  form, 


H(J,0;z)  =  H  (J)  +  Y,  Qj(J)  I  dk  a  (k  )  cos  [n0  ±  S(k  )z]  ,  (35) 

0  Yr  n=0  n  J  s  s  s  s 

where  &(ks)  is  given  by, 


&(k  )  =  &  +  (k  -  k  )  *-= 

v  s'  so  '  s  so'  k 


v(k  )  -  v(k  ) 


so 


(36) 


and  5so  =  ^kw/kr^kso  "  ^r ^  in  the  sPirit 
parenthesis  in  the  right-hand  side  of  (36) 

resulting  from  the  dispersive  effects  in  the 


of  Eq.  (10).  The  last 

2 

is  of  order  (w  /»  ) 
p  r7 

sideband  spectrum.  The 


27 


finite  ky  provides  phase  slippage  among  the  ponderomotive  phases  of 

various  wavenumbers  kg,  a  necessary  condition  for  the  validity  of  the 

quasilinear  theory,  even  when  the  dispersive  effects  of  the  last  term  in 

2  2 

(36)  are  negligible,  i.e.,  y  (w  /w  )  «  Dk  /k  . 

z  p  r  s  s 

The  resonant  condition  between  a  sideband  wavelength  kg  and  a 

given  harmonic  n  now  reads  +nub(Jn)  =  6(ks).  For  each  harmonic  n  there 

exists  a  wide  band  of  resonant  orbits  centered  around  J  and  of  width 

no 

DJ  defined  by, 
n 

DJn  .  (dV<U>Jno'1  (kv/kr>  Dks’  in“b<J„„>  *  5so-  <37> 

One  condition  for  the  applicability  of  the  quasilinear  theory  is 
that  the  phase  mixing  due  to  Dkg  occurs  much  faster  than  the  bounce 
period  around  a  secondary  island  in  phase  space.  In  this  way,  electrons, 
that  otherwise  would  execute  periodic  orbits  around  some  fixed  point, 
lose  coherence  sufficiently  fast  to  allow  random  motion  of  the  Fokker- 
Planck  type.  Taking  the  decorrelation  length  for  the  phase  1^  ~ 
2n/D6(ks)  and  applying  Eq.  (18a)  for  the  bounce  period  around  the  nth 
harmonic  we  obtain 


1/2 

I  <1 


i  n ) 


(  Aa  a 
w  r 


1/2 


U  +  a 


w 


(38) 


Inequality  (38)  guarantees  the  diffusive  behavior  within  the  band  DJ 

n 

around  JnQ  given  by  (37).  Large  scale  diffusion,  permitting  transport  of 
deeply  trapped  particles  across  the  separatrix  Jg  of  the  original 
primary  island,  requires  that  different  stochastic  bands  touch  each 
other,  DJn  +  DJn  +  ^  >  &Jn,  or,  using  a  similar  approach  as  in  Eqs.  (19) 
and  (20), 
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Q 

f  nQ  *1 

(  Aa  a 

_ s  .  _ 

n 

w  r 

k  0 

s  8n 

.  /~2  . 

.1  +  a  1 . 

rj  ' 

1/2 


(39) 


The  right-hand  sides  of  Eqs.  (38)  and  (39)  are  of  the  same  order  as  the 
upshifted  synchrotron  period  for  the  main  bucket.  Thus  both  conditions 
are  satisfied  when, 


»  2 


C 


v°> 


(40) 


Note  that  (40)  is  the  same  as  the  condition  (26)  that  defines  the  wide 
spectrum,  obtained  in  the  previous  subsection  using  different  arguments. 
Then  the  evolution  of  the  initial  distribution  fQ(J)  is  globally 
described  by  a  diffusion  equation, 


3f 

3z 


=  It  Vj)  If 


(41) 


20  21 

Applying  the  standard  procedures  of  the  quasilinear  theory  ’  (see 
appendix  C)  and  taking  the  limit  of  small  growth  rate  for  the  sidebands, 
Im(ks)/ks  <<  1,  we  obtain, 


k.  k  3.  oo 

Dq(J)  =  ~\T  "T  E  "2  |Qn±CJ)  I2  f  dk  V  (k  )  6(k  -  k  ), 
q  4  Kw  Y  n=0  n  J  s  s  b  s  n 

r 

k 

k„"  kr  1  IT  n  Kb<J)  ' 

W 


(42) 


According  to  the  condition  (40)  for  the  validity  of  the  quasilinear 
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theory,  the  wave  package  has  a  wide  spectrum  Dkg  ~  N(kr/kv)a)jj(0)  with  N 

large.  The  wave  components  ag(kn)  fall  off  slowly  for  n  <  N  while  the 

Fourier  coefficients  Qn  (J)  decay  rapidly  with  n.  Then  we  may  factor  out 

2 

the  average  spectral  power  density  V  =  (1/Dk  )J  dk  W  ( k  )  =  a  /Dk  in 

s  s  s  s  s  s  s 

Eq.  (34)  getting, 


VJ>  %-r  VsE,  "2IC<j>I2. 

H  4  w  y  n=0 
r 


The  infinite  sum  in  the  right-hand  side  of  (43)  is  computed  in  Appendix 
C.  The  summation  technique  does  not  require  the  knowledge  of  the 
individual  coefficients  Qn(J)  and  the  result  depends  only  on  the 
quantities  J  and  <*^(J)  for  the  unperturbed  Hamiltonian  Hq(J).  We  then 
obtain  the  diffusion  coefficient  in  closed  form, 


Dq(J)  = 


2  2 

k  a  y  V  k  J 
r  w  z  s  w 


VVJ> 


We  note  in  passing  that  the  method  used  to  obtain  expression  (44)  for 
Dq(J)  is  quite  general  and  valid  for  any  integrable  dynamical  system 
Hq(J)  that  is  subject  to  an  external  perturbation.  In  particular,  it 
should  be  applicable  to  a  variety  of  RF  heating  methods  in  fusion 
plasmas,  commonly  involving  a  strong,  narrow-band  pump  wave  embedded  in 
a  wide,  parametrically  excited,  fluctuation  spectrum. 

Using  the  expressions  (A3)  and  (A7)  for  the  action  J  and  the 
synchrotron  frequency  (^(J)  ve  find  that  the  diffusion  coefficient  goes 
to  zero  at  the  centre  of  the  primary  island  J  =  0,  has  a  logarithmic 
singularity  at  the  separatrix  J  =  Jg  and  falls  off  away  from  it.  In 
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normalized  units,  with  the  wiggler  wavelength  as  the  unit  length  and 
the  action  Jg  at  the  separatrix  as  the  unit  action,  we  have, 


D 

-q 


pq<J> 


(45) 


Choosing  the  value  J  =  J  /2  inside  D  (J)  we  obtain  an  estimate  for  the 

s  Q 

effective  diffusion  coefficient  associated  with  the  uniformly  filled 
distribution, 


D 

-q 


it  g  a  a  a 

&  y  r  _ s _ 

4  A  (1+a  2)  Dk  a  2 
w  s  r 


(46) 


where  g  is  the  same  scaling  factor  as  in  Eq.  (32). 

Note  that  both  expressions  (32)  and  (46),  corresponding  to  the 
two  different  spectral  types,  are  independent  of  yr*  Thus,  for  the  same 
wiggler  parameters  and  total  sideband  power,  the  detrapping  distance  in 
wiggler  periods  is  independent  of  the  electron  beam  energy.  The 
dependence  of  the  effective  diffusion  rate  on  the  wiggler  tapering 
enters  through  the  form  factor  g( q/  )  =  J  (n)/J(*  ).  As  the  rate  of 
tapering  increases  and  shifts  further  from  n,  g  inreases1  and 
accelerates  the  effective  diffusion  rate.  This  happens  because  the 
trapped  area  in  phase  space,  parametrized  by  Jg,  is  shrinking  as  the 
tapering  progresses,  while  the  sideband  induced  excursions  remain  the 
same,  depending  mainly  on  the  sideband  strength  and  spectrum.  This 
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shortens  the  average  detrapping  time  for  a  particle.  The  diffusion  by 

2 

broad  discrete  spectrum,  Eq.  (32)  scales  as  g  ,  while  the  quasilinear 
diffusion,  Eq.  (46),  scales  as  g.  Thus  the  former  is  affected  more  by 
tapering  than  the  latter. 
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VII.  NUMERICAL  RESULTS 


The  numerically  computed  diffusion  coefficient  D  and  the 

diffusion  length  in  wiggler  periods  =  1/D  are  plotted  against  the 

2  2 

sideband  to  carrier  power  ratio  P  =  E  a  («  )/a  =  V  /W  in  Figs.  7  and 

ns  n  r  sr 

8  respectively  for  the  three  different  types  of  spectra.  We  have 

integrated  numerically  the  equations  of  motion  for  400  particles  of  a 

uniform  initial  distribution  inside  the  bucket.  The  field  intensities 

remained  constant  at  a  =  5x10  a  =2  and  y  =  25.  All  the  numerical 

r  v  r 

results  in  this  paper  correspond  to  a  tapered  wiggler  with  ^  =  7 n/6 .  A 
clear  separation  in  the  diffusion  rates  is  observed  among  the  various 
spectral  types.  The  narrow  frequency  results  (triangles)  were  obtained 
using  the  Hamiltonian  (8)  with  a  single  sideband  frequency  »  /w  = 

1.016.  The  results  for  a  broad  discrete  spectrum  (circles)  were  obtained 
using  (23)  with  A  =  20,  «  /«  =  1.016  and  v  =  0.5  S  .  The  continuous 

S  L  S 

spectrum  (squares)  was  modeled  by  A  =  100,  v  =  0.05  5  .  The  solid  lines, 

corresponding  to  the  theoretical  results  of  Eqs.  (32)  and  (46),  are  in 

good  agreement  with  the  numerics.  Theoretical  predictions  for  the  single 

frequency  case  were  not  made.  We  stress,  however,  the  difference  between 

single  frequency  results  and  quasilinear  theory  in  this  case.  The 

22 

agreement  that  has  been  observed  in  some  other  cases  is  not  generic 

but  particular  to  certain  systems. 

Figure  9  illustrates  the  difference  of  the  electron  response  to 

different  sideband  spectra.  The  plots  on  the  left  side  are  typical 

orbits  J(z)  for  selected  particles  along  the  wiggler.  The  trajectories 

in  all  plots  are  generated  by  the  same  initial  conditions  for  the 

electrons  and  the  same  FEL  parameters  aT  ,  a„  and  k  ,  as  well  as  the  same 

v  r  w 

2 

mean  square  sideband  power  <ag  >.  The  spectral  parameters  A  and  v, 
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however,  are  different  so  that  each  of  the  figures  (a)  to  (c) 
corresponds  to  one  of  the  three  spectral  types  defined  earlier.  The 
dashed  line  marks  the  position  of  the  unperturbed  separatrix  Jg.  The 
corresponding  distribution  functions  f(J,z)  at  the  beginning,  z  =  0, 
halfway  inside,  z  =  50Xw,  and  at  the  end,  z  =  100Xw  of  the  wiggler  are 
plotted  in  the  right-hand  side  of  Figs.  9a-9c  respectively. 

In  Fig.  10  we  plot  the  diffusion  coefficient  for  a  uniformly 
filled  bucket  as  a  function  of  the  energy  vr»  fixing  the  wiggler 
parameters.  It  is  clear  that  the  diffusion  rate  (measured  again  in 
number  of  wiggler  periods)  is  independent  of  the  beam  energy,  provided 
the  synchrotron  frequency  stays  in  the  same  parameter  regime. 

Once  the  diffusion  coefficients  are  known,  some  estimate  can  be 
made  of  the  related  reduction  in  efficiency  over  the  wiggler  length.  The 
number  of  trapped  particles  at  any  point  z  is  given  by 

n.  =  \S  dJ  f(J,z).  Using  the  diffusion  Eq.  (41)  with  D(J=0)  =  0 
n  J0 

one  obtains  the  rate  of  change  in  the  number  of  trapped  particles, 
dn.  3f(J  ,z) 

d z  =  Vz>  D(V  —r—  •  (47> 

dZ 

The  leakage  rate  for  trapped  particles  changes  along  z  as  the  slope  of 
the  distribution  f(J,z)  changes.  To  estimate  the  average  leakage  rate  we 
assume  that  f(J,z)  remains  Gaussian  in  J  with  an  average  width  equal  to 
the  separatrix  action  Js.  We  estimate  from  (47)  the  e-folding  length 
=  -  n^  ^(dn^/dz)  for  the  number  n^  of  trapped  particles, 


Ld  =  D<V  Js 


-2 
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(48) 


Assuming  that  nb(z)  =  nb<0)  exp(-z/Ld),  the  number  of  detrapped 

particles  between  z  and  z  +  Az  is  An(z)  =  nb(0)Ld~*  exp(-z/Ld)  A z.  These 

particles  gave  up  an  amount  of  energy  AE(z)  =  [rr(0)  -  rr(z))  An(z)  as 

radiation.  Integrating  AE(z)  over  the  wiggler  length  for  a  linearly 

tapered  wiggler  Yr(z)  =  rr(0)  -  zAy/Lv,  we  find  the  total  energy 

extracted  from  the  particles  that  were  detrapped  at  some  point  inside 

the  wiggler.  Adding  the  contribution  (y(0)  -  y_(L)]  n,  (L  )  from  the 

r  r  v  D  w 

particles  that  remained  trapped  throughout  the  wiggler  length,  we  come 
up  with, 


h  = 


1 


exp  (- 


(49) 


where  =  Ay/Yj.(0)  is  the  efficiency  without  induced  diffusion.  The 
loss  of  amplification  will,  in  general,  be  distributed  among  all  the 
radiation  modes  and  (49)  reflects  the  total  power  loss  in  all 
frequencies.  The  extraction  efficiency  t)  for  a  linearly  tapered  wiggler 
is  plotted  in  Fig.  11  versus  the  sideband  to  carrier  power  ratio  P, 
obtaining  the  corresponding  value  for  from  the  results  in  Fig.  8. 
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VIII.  CONCLUSION 


The  diffusion  in  phase  space  caused  by  sideband  excitation 

during  FEL  operation  was  studied.  It  was  shown  that  the  characteristic 

rates  for  this  process  depend  on  the  structure  of  the  sideband  spectrum, 

falling  into  one  of  the  following  general  categories:  narrow,  wide 

discrete  or  wide  continuous  spectrum.  In  all  cases,  the  diffusion 

coefficient  was  found  proportional  to  the  ratio  of  the  total  power  in 

the  sidebands  to  the  power  in  the  main  FEL  signal  D  *  C  Vs/V  .  The 

coefficient  C,  however,  is  connected  to  the  spectral  type  under 

consideration.  From  Eqs.  (32)  and  (46)  we  see  that,  apart  from  numerical 

factors  of  order  unity,  C  scales  as  (a  a  )^  with  1  =  1/2  for  a  discrete 

v  r 

and  1=1  for  a  continuous  spectrum.  Therefore,  given  the  typical  FEL 

_3 

values  of  a  <10  and  a  <  10  an  order  of  magnitude  reduction  in 
w  r  - 

diffusion  occurs  in  the  transition  from  a  discrete  to  a  continuous 
sideband  spectrum.  It  was  also  observed  numerically  that  the  highest 
diffusion  rate  occurs  when  all  the  sideband  power  is  in  a  single 
frequency.  In  this  case,  however,  a  portion  of  the  particles  will  remain 
trapped  for  arbitrary  long  wigglers  if  the  sideband  amplitude  is  below 
the  threshold  for  extensive  stochasticity.  The  stochasticity  threshold 
is  progressively  reduced  as  the  sideband  power  is  distributed  into  an 
increasing  number  of  frequencies.  Yet  the  rate  of  diffusion  also  slows 
down  with  increasing  spectral  width  and  decreasing  mode  separation. 

Thus,  the  minimum  reduction  in  the  FEL  energy  extraction  efficiency  will 
occur  for  continuous  sideband  spectra.  Although  control  of  the  sideband 
structure  does  not  seem  plausible,  experiments  show  that  a  wide  spectrum 
is  naturally  excited  during  FEL  operation.  This  would  allow  enough  power 
build-up  before  serious  deterioration  in  efficiency,  due  to  detrapping, 


36 


to  occur.  The  diffusion  length,  measured  in  wiggler  periods,  is 
independent  of  the  beam  energy  y  under  the  same  wiggler  parameters,  for 
all  the  spectral  types.  Our  results  have  been  obtained  for  radiation 
fields  of  constant  amplitude.  Inclusion  of  the  time  evolution  for  both 
the  carrier  and  the  sidebands  will  modify  the  detrapping  rates  by 
changing  both  the  diffusion  rate  as  well  as  the  size  of  the  separatrix. 
This  subject  is  left  for  future  investigation. 


ACKNOWLEDGEMENT 


This  work  was  supported  by  SDIO  and  managed  by  SDC. 


APPENDIX  A:  TRANSFORMATIONS  IN  ACTION-ANGLE  VARIABLES 


The  relations  between  y,  ♦  and  the  action-angle  variables  J,  © 
are  given  in  closed  forms  in  case  of  an  untapered  viggler.  Starting  from 
the  general  expression  (11)  and  using  (8)  we  have, 


9  = 


-  mx 

r  v  h 

1/2 

d<J» 

-  V(*) 

» 

*mn 

rrKb(J) 

l  k 
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d<J/' 
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J 

r  y  h 

r  c 

a  a 

»  w  r 

2k 

w  “ 

„  k 
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kw 

,-1/2 


V(Y) 


where 


V(*J»)  »  ^  cos<J*  +  <J»  sinij/^  j  . 


(Al) 


(A2) 


Using  H  =  -  (a  a  /y  )  V(<|>  )  we  obtain  the  action  at  the  separatrix, 

S  W  S  L  L 


a_,l/2  1 


r  mx 


y  n  J 

r 


1/2 


d*  (  V(*)  -  V(*r)  ] 


mn 


In  case  of  an  untapered  wiggler  \|»r  =  0  Eqs.  (A1)-(A2)  yield, 


Jg  (e2<\)  -  (1-X2)E1(X))  ,  X2  <  1 


J  = 


2JsX  E2  IxJ  ’ 


2 

XZ  >  1 


(A3) 
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Aawackt 


,1/2 


4k 


w  J 


(A4) 
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X2  <  1, 


Xsn^E1(X)e)  , 
sn  (n  El(x)9)  - 


(A5) 


where  E^  and  E£  are  the  complete  elliptic  integrals  of  the  first  and  second 
kind,  sn  is  the  Jacobi  elliptic  sine  function,  and 


X*  = 


Y  H 
r  0 

2a  a 
w  r 


+  — 


(A6) 


is  the  trapping  parameter  (  XZ  <  1  for  trapped  particles).  Using  ^(J) 
=  (3J/3Hq)_*  and  (Al)  we  determine  the  bounce  (synchrotron)  frequency 


Kb(J)  - 


Kb(0)  2E*(X)  ’ 

itX _ 


Kb(0)  Ex(l/X) 


X2  <  1 


X2  >  1 


(A7) 


where 


v°> 


1_ 

Y. 


(7  a  a  k  k  ] 
(2  v  r  w  rj 


1/2 


(A8) 


2 

is  the  bounce  frequency  at  the  bottom  of  the  well.  J  and  X  are  mutually 
related  through  (A3)  -  (A6)  and  they  uniquely  label  the  trajectories. 

The  Fourier  coefficients  of  the  expansion  (12)  can  also  be 
expressed  in  closed  form.  They  are  computed  by  integration  in  the 
complex  plane  around  the  singularities,  utilizing  the  double  periodicity 
properties  of  the  Jacobi  elliptic  functions  to  obtain, 


Qn  = 


-<±1)‘ 


nn 


n 

2 


E^(X) 


l-(-q) 


q  =  exp 


rnE^xn 


*X>  J 


2 

XZ  <  1, 


(A9) 
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+IO 


2  2 

n±  nn  X  n 
Qn  =  -  IT -  « 


e;(i/x) 


i  i 

..  2n  ±  . 
1-q  1+q 


2n 


q  =  exp 


|'JlE1(l/X)'| 


[  Ejd/XJJ 


,  X*  >  1 


for  n  t  0  and, 


=  2 


E2(X) 

E^X) 


-  1  XZ  <  1  ,  Qq  =  1  -  2X2 


E^l/X)  J 


where  e'(X2)  a  E1(l  -  X2), 
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APPENDIX  B  :  PHASE  AVERAGING  OVER  CONSTANT  J 


The  phase  average  <  ! v ( >  ~  V(^)  j  >  over  constant  J  is  given 
for  i pr  =  0,  by, 


r2n 


>  =  — 
2n  J 


d9.  |  cos^mx(J)  -  cos^CJ,^)  | 


4 

2n  J 


vj/ 

/mx 


d*i  fc)  I  C0S*mx(J)  -  cos*i  I 


(Bl) 


Substituting  d8  /d<J/.  from  (A2)  and  using  cos  <f»  =  y  H(J)/a  a  one 

n  x  mx  r  v  r 

obtains, 


.>  = 
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(  MX  >  2 

— — j-  -  (  1  -  XZ) 
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(B2) 


where  X^  was  defined  in  (A6). 


APPENDIX  C  :  COMPUTATION  OF  THE  QUASI LINEAR  DIFFUSION  COEFFICIENT 


Ve  consider  the  evolutit i  of  the  electron  distribution  in  the 
presence  of  a  sideband  wave  package. 


as<z’‘>  ■  n  J 


dk  a  (k  ) 
s  s  s 


ikgz  -  i»(ks)t 


(Cl) 


of  width  Dkg.  The  interaction  Hamiltonian  in  action-angle  variables, 
derived  in  Sec.  IV. (b)  is, 


a  ®  - 

H(J, 0;z)  =  H  ( J)  +  -*•  £  Q;(J>  dk  a  (k  )  cos  [n9  ±  S(k  )z]  . 

°  Yr  m=0  n  J  s  s  s  s 

(C2) 

The  evolution  of  the  distribution  function  f(J,6;z),  under  the 
Hamiltonian  flow, 


d0  _  3H  dJ  _  9H 

dz  "  3J  ’  dz  "  39  ’ 


is  given  by, 


(C3) 


3f  de  3f  dJ  3f  _ 
3z  +  dz  30  +  dz  3J 


(C4) 


Ue  separate  the  distribution  f(J,0;z)  into  a  slowly  varying  part  fQ(J;z) 
=  <  f  >  and  a  fluctuating  part  Sf(J,0;z)  =  f  -  <  f  >.  The  averaging 
operator  is  defined  by, 


<  f  > 
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It  is  implied  in  the  above  definition  that  the  characteristic  length  for 
f o( J ; z )  is  longer  than  the  synchrotron  length  L  =  2 n/ .  We  then  obtain 
from  (C4), 


af  deaf  dJ  3f  _dj!£o_a_rdj: 
az  +  dz  ae  +  dz  aj  ~  dz  aj  aj  [  dz  r 


]  ,  (C5) 
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o 


(C6) 


Using  Eqs .  (C2)  and  (C3)  inside  Eq.  (C5)  and  ignoring  the  last  bracketed 
term  in  the  right-hand  side  we  obtain, 


f  =  -  —  E 


2 Yr  n=lj 


n  Q  ~(J)  a  (k  )  3f  .  r  Q 
,,  n  sv  s'  o  i[n0  +  &fk  )1 

dk - e  1  s'  J  +  cc  . 

s  n  Kb(J)  ±  S(ks)  3J  (C7) 


Substituting  (C7)  in  (C6)  we  have, 
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n  Qn*(J)  as(ks)  3fo 
n  Kb(J)  ±  8(kg)  3J 


ei[n0  t  S(kg)J 


9 


where  again  5(kg)  =  (kw/kr) { (kg-kr )  +  i e] ,  e  =  Im(ks).  Integration  of  the 
right-hand  side  of  (C8)  over  0  yields, 


IQ2  as*(qs)  as(ks) 
n  Kb(J)  ±  S(ks) 


&(k  )Jz 

+  cc 
(C9) 
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Spatial  integration  yields 


APPENDIX  D:  SUMMATION  OF  FOURIER  COEFFICIENTS 


We  present  a  general  technique  of  computing  sums  of  the  form, 


E  n2  [  |Qn+(J)|2  +  |Qn“(J)|2  ]  •  (Dl) 

n=0  L  n  n  J 

The  quantities  Qn-(J)  are  the  Fourier  coefficients  from  the 
decomposition  of  the  phase  exp[  i<|>( J, 0)  ]  of  the  perturbation  into 
harmonics  of  the  angle  variable  0  for  the  unperturbed  system.  The 
knowledge  of  the  individual  Qn~(J)  *E  not  required  in  the  computation. 
The  technique  should  be  applicable  to  a  wide  class  of  integrable 
systems  experiencing  a  periodic  perturbation  with  only  minor 
modifications.  In  our  case  Qn~(J)  are  defined  by, 


cos  [«KJ,0)  +  &sz] 


00 

Y  Q+( J)  cos(n0  +  6  z)  +  Q~(J)  cos(n0  -  6  z). 
n=0  n  s  n  s 

(D2) 


Closed  forms  for  Qn-(J),  obtained  in  Ref.  2  for  the  case  of  an 
untapered  wiggler,  appear  in  Appendix  A. 

For  untrapped  particles  we  have, 


Q2m'  =  Q2m  ’ 


°2m+l  =  4  Q2m+1 ' 


Setting  &s  =  0  in  (D2)  and  differentiating  in  0,  we  obtain, 


gg  cos  <p 


d<j> 

d0 


sin  {//  = 


(2m)  2Q2m  sin  2m0  , 


(D3a) 
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sin  I9  '  fl  =  Jg  cos  <1  =  E  (2m+l)  2Q2m+1  cos[ (2m+l)d) . 

"'°  (D3b) 


Squaring  the  right-hand  sides  of  (D3a),  (D3b),  adding  them  and 
integrating  over  9  we  obtain, 


E"2  [|Q„*|2  *  |0„-|2]  .  2  E"2|Qn|2-k 


^2n 


m}2  (cosZ'<'  +  sin2^)  dO 


(D4) 


Applying  the  same  procedure  to  untrapped  particles  we  obtain, 


30  C0S  *  =  -  d0 


,  °° 

sin  <J/  =  -  £  n  (Qn+  +  Qn_)  sin  n0 


m=l 


(D5a) 


CD 

|q  sin  [0  -  f]  =  cos  q  =  -  £  n  (Qn+-  Qn  ~)  cos  n0  ,  (D5b) 


2 J  ~  d0 


m=0 


n  n 


and,  after  squaring,  adding  (D5a)  and  (D5b)  and  integrating  over  0  , 
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Thus  in  both  cases, 


nE  "2  (|Q„*|2  *  |Qn-|2) 


l_ 

2n 


d0 


(D7) 


Using  the  definition  Eq.  (11)  for  0(t),  and  Eq.  (8)  for  the  unperturbed 
Hamiltonian  (ag  =  0  )  in  the  right-hand  side  of  (D7)  we  have, 
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The  last  integral  in  (D8)  is  by  the  definition  (11)  the  action  J  for 
the  unperturbed  Hamiltonian,  yielding  the  final  result, 
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Figure  6.  The  fraction  fd  of  detrapped  particles  after  100  viggler 
periods  for  an  initially  uniformly  filled  bucket.  Results 


for  various  sideband  frequency  ratios  m  =  c^/co^  are 
plotted  against  the  relative  sideband  amplitude  q  =  ag/ar 
for  &  =  2,  a„  =  5xl0-^  and  *  =  7n/6,  y  =  25. 


Figure  7.  The  normalized  diffusion  coefficient  D  as  a  function  of 

the  sideband  power  ratio  P  =  for  ay  =  2,  ar  =  5xl0-^ 

and  =  7n/6,  rr  =  25.  Squares  correspond  to  a  continuous 

type  of  spectrum,  peaked  at  u  /»  =  1.016  with  A  =  100  and 

v  =  0.1.  Dots  correspond  to  a  wide  discrete  peaked  at 

ws/«r  =  1.024  with  A  =  20  and  v  =  0.5.  Triangles 

correspond  to  a  single  frequency  spectrum  with  »  /«  = 

s  r 

1.016.  The  upper  and  lower  solid  lines  correspond  to  the 
theoretical  results  from  Eqs.  (36)  and  (42)  respectively. 


Figure  9.  Particle  response  under  different  sideband  spectra.  On  the 

left-hand  side  we  plot  the  action  J  against  the  number  of 

wiggler  periods  Ny  for  selected  particles.  On  the  right- 

hand  side  we  plot  the  corresponding  distribution  function 

f(J)  at  Ny  =  0,  50,  and  100.  In  all  cases  the  total 

sideband  power  ratio  P  =  W  /W  =  0.36  and  a  =  2,  a  = 

s  r  v  r 

5xl0~5,  vr  =  25.  (a)  corresponds  to  a  wide  continuous 

sideband  spectrum  (b)  corresponds  to  a  vide  discrete 
spectrum  and  (c)  to  a  single  frequency  «  /«  *  1.016. 

S  TO 
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Plots  of  the  normalized  diffusion  coefficient  D  against 

the  electron  energy  yr*  Parameters  are  ay  =  2,  ar  =  2  x 

10'4,  a  =  7.5  x  10~5  and  <o  =  1.024.  Squares  correspond 
s  s 

to  a  continuous,  dots  to  a  wide  discrete  and  triangles  to 
n  single  frequency  spectrum. 


u/u 


1.2 


Lw=  100  X 


w 


P 

Figure  11.  Deterioration  of  the  extraction  efficiency  for  various 
spectral  types.  Plotted  is  the  ratio  for  a  viggler 

length  of  100  periods  as  a  function  of  the  sideband  power 
ratio  P.  The  numerical  results  in  Fig.  8  were  used  in  Eq. 
(49)  to  create  this  plot.  The  symbolism  and  the  parameters 
are  the  same  as  in  Fig.  8. 
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